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Abstract 

Two sets of vectors, covariant and orthogonal Lyapunov vectors (CLVs/OLVs), are currently 
used to characterize the linear stability of chaotic systems. A comparison is made to show their 
similarity and difference, especially with respect to the influence on hydrodynamic Lyapunov modes 
(HLMs). Our numerical simulations show that in both Hamiltonian and dissipative systems HLMs 
formerly detected via OLVs survive if CLVs are used instead. Moreover the previous classification of 
two universality classes works for CLVs as well, i.e. the dispersion relation is linear for Hamiltonian 
systems and quadratic for dissipative systems respectively. The significance of HLMs changes in 
different ways for Hamiltonian and dissipative systems with the replacement of OLVs by CLVs. 
For general dissipative systems with nonhyperbolic dynamics the long wave length structure in 
Lyapunov vectors corresponding to near-zero Lyapunov exponents is strongly reduced if CLVs are 
used instead, whereas for highly hyperbolic dissipative systems the significance of HLMs is nearly 
identical for CLVs and OLVs. In contrast the HLM significance of Hamiltonian systems is always 
comparable for CLVs and OLVs irrespective of hyperbolicity. We also find that in Hamiltonian 
systems different symmetry relations between conjugate pairs are observed for CLVs and OLVs. 
Especially, CLVs in a conjugate pair are statistically indistinguishable in consequence of the micro- 
reversibility of Hamiltonian systems. Transformation properties of Lyapunov exponents, CLVs and 
hyperbolicity under changes of coordinate are discussed in appendices. 
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I. INTRODUCTION 



Chaos means a sensitive dependence on initial conditions. This intrinsic randomness 



_ treatment of them feasible, which is 
l|. Besides, chaos plays an important 



of the fully deterministic systems makes a statistica 
essential for the foundations of statistical mechanics 
role in a plenty of phenomena which is of relevance to our daily life, for instance the weather 
forecasting 

To characterize the chaoticity of dynamical systems, Lyapunov exponents and vectors are 
mostly used. One important recent finding of Lyapunov analysis is that for systems with 
continuous symmetries Lyapunov vectors corresponding to near-zero Lyapunov exponents 
have long wave-length structures, named hydrodynamic Lyapunov modes (HLMs) 3|. This 
provides a new possibility to connect the reduced description of a many-body system to the 
microscopic information of its detailed dynamics. Further investigations showed that HLMs 
exist in a large number of systems and they have some universal features irrespective of 
the details of their dynamics js]. One should mention that localization of Lyapunov vectors 
corresponding to the largest Lyapunov exponents was also intensively studied jjj. 

Lyapunov analysis was conventionally undertaken via the so-called Benettin algorithm 



111 ] . where Lyapunov vectors are calculated as the set of orthogonal vectors right after 



reorthogonalization of offset vectors. Recently, the application of another set of vectors called 
covariant Lyapunov vectors (CLVs) was made feasible via an efficient algorithm proposed by 
Ginelli et al. [U] . CLVs have been shown suitable for the characterization of hyperbolicity 
of high dimensional systems since they are expected to span the local stable and unstable 
subspaces of the investigated systems. In view of the obvious difference, it becomes necessary 
to study the relation between CLVs and the conventionally used Lyapunov vectors. We 
denote the latter as orthogonal Lyapunov vectors (OLVs) in order to distinguish them from 
CLVs. 

We first recall in Sec. II the definition and numerical calculation of both sets of vectors. 
The model system of coupled map lattices (CMLs) is introduced in Sec. III. Through inten- 
sive numerical simulations the following questions are addressed in the remaining sections: 
(i) will HLMs survive if CLVs are used instead of OLVs (Sec. IV), (ii) are HLMs from CLVs 
as significant as those from OLVs, (iii) what is the implication of the Hamiltonian structure 
to the relation between conjugate pair of CLVs, and what is the implication for the relation 
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between coordinate and momentum parts of CLVs? 



II. DEFINITION AND CALCULATION ALGORITHM FOR CLVS AND OLVS 



Recall that in the seminal work 10) about the multiplicative ergodic theorem Oseledec 
proved that the limit H = lim^ +00 [M T (t, 0) • M(t, 0)]» exists for almost every initial point 
of a nonlinear dynamical system, where M(t, 0) is the fundamental matrix governing the 
time evolution of perturbations 8X(t) in tangent space as SX(t) = M(t, 0) • SX(0). The set 
of Lyapunov exponents are defined as X^ = ln/x^ a \ where /z^ are the eigenvalues of the 
matrix S, i.e. S • = fx^g^ a \ 

In practice the Lyapunov exponents and OLVs are calculated via the so-called standard 
method invented by Benettin and Shimada et al. which was used in most studies of 
HLMs Here the time evolution of a set of offset vectors in tangent space is moni- 

tored by integrating the linearized equation. And the offset vectors are reorthonormalized 
periodically. The time averaged values of the logarithms of the renormalization factors 
are the Lyapunov exponents and the set of offset vectors right after the reorthonor- 
malization are the Lyapunov vectors. The relation between the Oseledec eigenvectors g^ 
and the Lyapunov vectors obtained via the standard method is subtle. It is proved that 
the Lyapunov vectors f w obtained via the standard method converge exponentially to the 
Oseledec eigenvectors for the inverse-time dynamics of the original system 
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13J. In 



other words, as t — > +oo, there is ~ g^ a \ where g^ are eigenvectors of the matrix 



lim^ +00 [M(t, 0)' ■ M(t, 0)]* as well as its inverse 5 * = lim t ^ +00 [M(t, 0) ■ M(t, 0) T ]s 



and M(U0) = [M(t, 0)] 1 is the fundamental matrix of the inverse-time dynamics. See 



Rcf. 
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13] for the details. 



Exactly in the same work 10| Oseledec proved also that, for almost all initial conditions 
x, there is a splitting of the tangent space TM(x) 



TM(x) = E\x)@E 2 (x) 



E s (x), 



(1) 



and there exist real numbers Xi(x) > X^ix) > ■ ■ ■ > X s (x) such that 



lim -]ndDf n \ E id = Xi(x), 

n->±oo n 



(2) 



where Df is the derivative governing the tangent space dynamics. The set of numbers Xi(x) 
with degeneracy rrii = dimE^x) composes the Lyapunov spectrum and the decomposition 



stated in Eq.([T]) is called the Oseledec splitting. The spanning vectors e( a ) of the Oseledec 
subspace E l {x) are the CLVs. 

In contrast to the popularity of OLVs, the use of CLVs was made feasible only recently 
owing to an efficient algorithm proposed by Ginelli et al. jl^j]. The new algorithm relies 
on the information obtained via the standard method of Benettin. One additional integra- 
tion of the inverse-time dynamics is performed in order to get CLVs and the corresponding 
fluctuating finite-time Lyapunov exponents. The basic idea is that an arbitrary offset vec- 
tor will approach asymptotically the most unstable direction corresponding to the largest 
Lyapunov exponent. It is known that the covariant fc-dimensional subspace spanned by the 
first k CLVs is spanned by the first k OLVs as well. An arbitrary offset vector confined to 
this subspace will approach asymptotically the fc-th CLV if the inverse-time tangent space 
dynamics is applied. To this aim one needs the effective tangent space dynamics confined 
in the /c-dimensional covariant subspace. A representation of this effective dynamics in the 
coordinate space of OLVs is given by the R-matrix produced by the reorthonormalization 
steps of the standard method. Detailed formulas can be found in Ref. [ijj]. 



151 ] . It is demanded to compare 



We mention that a different algorithm was used in Ref. 
the efficiency of the two algorithms. 

To characterize Lyapunov vectors of extended systems quantitatively, we introduced in [6j 
a dynamical variable called LV fluctuation density in the spirit of generalized hydrodynamics, 

L 

U^(r,t) = Y,^ { t a)l -5(r-r l (t)), (3) 
1=1 

where 5(x) is Dirac's delta function, ri(t) = I ■ a is the position coordinate of the l-th 
element taken here as ri(t) = I ■ a, and {8u[ } is the coordinate or momentum part of the 
a-th Lyapunov vector at the discrete time t. For simplicity, we set a = 1 in the following 
discussion. The spatial structure of LVs is characterized by the static LV structure factor 
defined as 

(k) = [ (U (a) (r, 0)W (Q) (0, 0))e- jk r dr, (4) 



which is just the spatial power spectrum of the LV fluctuation density. 

As shown in past studies jfl-sj], the quantity k max representing the wave-number of the 
dominant peak of the structure factor S(k) and S(k max ) can be used to characterize the 
significance of long wave length structure in Lyapunov vectors. 
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III. MODELS 



CMLs [16| were selected as the main focus of this study because they have, which is 
essential to HLMs, similar symmetries as many-particle systems but are relatively much 
simpler. 

The two classes of CMLs under investigation have the form 



v 



i 

t+i 



v\ + e[f{u\^-u\)-f{u\-u\- 1 )] (5a) 



< + i = < + 4+1 (5b) 

and 

<i = u\ + e[f(u l t +1 - u\) - f{u\ - rij" 1 )] (6) 

where f(z) is a nonlinear map, t is the discrete time index, / = {1, 2, • • • , iV} is the index of 
the lattice sites and N is the system size. Unless explicitly stated, we use periodic boundary 
conditions {u® = u^,uf +1 = u\} in the numerical simulations below. 

Two options of the local map are used, the sinusoidal map fc(z) = ^ sin(27rz) and the 
skewed tent map 

{z'/r for < z' < r, 

(7) 
{l-z')/(l-r) for r < z' < 1. 

with z' = z (mod 1). With the parameter being close to zero Eqs. ([5]) and ([6]) with the 
skewed tent map is highly hyperbolic whereas Eqs. ([5]) and ([6]) with the sinusoidal map is 
nonhyperbolic. Especially the Hamiltonian system Eq. ([5]) with the two options of the local 
map are similar to the well-studied cases of hard-core systems and soft-potential systems, 
respectively. Alternatively, tuning the parameter r of Eq. (E3) from to 0.5 leads to a smooth 
variation of the dynamics of Eq. ([5]) from hard-core-like to soft-potential like 17]. 

Obviously, both systems Eq. ([5]) and ([6]) are invariant under an arbitrary translation in 
u-direction. Such a symmetry is known to be responsible for the appearance of HLMs. 



IV. EXISTENCE OF HLMS IN CLVS 

We show in FigJTJthe contour plot of the static CLV structure factors S(k). For both 
dissipative and Hamiltonian systems, either with the skewed tent map or the sinusoidal 
map, a clear ridge structure can be seen in the regime (fc, A) ~ (0,0), which indicates 
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FIG. 1: Contour plot of the static CLV structure factors for (a,c) Hamiltonian and (b,d) dissipative 
coupled map lattices. The local dynamics used is the skewed tent map (a,b) and the sinusoidal 
map (c,d), respectively. Other parameters are e = 1.3 and r = 0.15. A ridge structure can be 
clearly seen in the small (A, k) regime which indicates the existence of HLMs. 



the existence of long wave-length structures in CLVs associated with near-zero Lyapunov 
exponents. These numerical results demonstrate that HLMs formerly detected via OLVs 
survive if CLVs are used instead. 
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FIG. 2: Dispersion relations X-k max obtained from CLVs and OLVs for (a) dissipative and (b) 
Hamiltonian systems, respectively. The local map is the skewed tent map with e = 1.3 and 
r = 0.15. Note the perfect agreement between data from CLVs and OLVs for the highly hyperbolic 
cases shown here. 



V. UNIVERSALITY OF DISPERSION RELATIONS 

In Ref.[3, B] we found that the X-k dispersion relation of HLMs can be classified into 
two universality classes with respect to the system dynamics. Dissipative systems have a 
quadratic X-k dispersion while Hamiltonian systems have a linear one. Now we see whether 
such a classification is still valid if CLVs are used instead. 

The cases with the skewed tent map as local dynamics are shown in FigfSJ As can be seen 
from the plot, the CLV dispersion relations for dissipative and Hamiltonian systems have 
the different asymptotic behavior. The former is of the asymptotic form A ~ ^ m while the 
latter is A ~ k max , as reported for OLVs [7]. Moreover, for the used parameter setting the 
dispersion curves for CLVs agree very well with those of OLVs . Cases with the sinusoidal 
map as local dynamics are shown in Figj3j For both dissipative and Hamiltonian systems 
CLV dispersions are converging to the expected asymptotic forms, even better than OLVs. 
Note also that, as shown in Figj3b, for Hamiltonian systems CLV dispersions for positive 
and negative Lyapunov exponents follow the same curve while OLV dispersions behave 
differently in the positive and negative Lyapunov exponent regimes. For the used system 
size only the positive Lyapunov exponent branch of OLV dispersion is close to the asymptotic 





• CLV; k>0 
« CLV; X<0 
OLV: \>0 

— OLV: X<0 

— y=4.3x 



FIG. 3: Similar to Fig. [2] but the local map is the sinusoidal map and e = 1.3. Note that for the 
nonhyperbolic cases shown in this figure CLVs and OLVs still have the same asymptotic behavior. 



form. Further discussion regarding these differences will be given in the following sections. 
Nevertheless, for the two representative cases the investigated CLV dispersions follow well 
the reported classification of the universality classes of HLMs [?], sj]. 



VI. SIGNIFICANCE OF HLMS 



To characterize the significance of long wave length structure in Lyapunov vectors, we use 
the measure S(k max ), which is the height of the dominant peak in the static LV structure 
factor S(k) (Eq. (jU). 

The dominant wave number k max and the significance measure S{k max ) are compared for 
CLVs and OLVs in FigfJ] for the cases with the skewed tent map as the local dynamics. For 
such highly hyperbolic systems both the position and the height of the dominant peak are 
nearly identical for CLVs and OLVs for either the dissipative system or the Hamiltonian sys- 
tem in the positive Lyapunov exponent regime. We postpone the discussion of the negative 
Lyapunov exponent part of the Hamiltonian system to the next section. This observation 
indicates that for the highly hyperbolic systems the significance of HLMs is not influenced 
if CLVs are used instead of OLVs. A similar comparison was made also for cases with the 
sinusoidal map as the local dynamics as shown in Fig. [5j For both, dissipative system 
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FIG. 4: Dominant wave number k max and the significance measure S(k max ) of CLVs and OLVs 
for (a,b) dissipative and (c,d) Hamiltonian systems respectively. The local map is the skewed tent 
map with e = 1.3 and r = 0.15. The long wave length structure is as significant in CLVs as in 
OLVs for the highly hyperbolic cases shown. 

and Hamiltonian system, clear discrepancies between CLVs and OLVS can be seen in k max 
and S{k max ). For the dissipative case, the height of the dominant peak S(k max ) for CLVs 
is much lower than for OLVs (see Figj5]), especially in the regime A ~ (a/N ~ 0.93), 
which means that for the strongly nonhyperbolic systems as shown here, the use of CLVs 
reduces the visibility of long wave length structure as compared to OLVs. In contrast, for 
the Hamiltonian case, the height of the dominant peak S(k max ) is comparable for CLVs and 
OLVs. Note, however, that the variation of S(k max ) for CLVs is symmetric with respect to 
the spectral center a/2N = 0.5 while it is asymmetric for OLVs (Figl5]i). 

To demonstrate further the influence of hyperbolicity on the significance difference be- 
tween CLVs and OLVs we tune the parameter r of the skewed tent map Eq.((7]). Results 
for dissipative cases and Hamiltonian cases are shown in Fig. [6] and [7] respectively. In con- 
sistence with our previous results in Ref. |l7| the weakening of hyperbolicity as increasing 




,0 " i zE^Oii — " i i i 

' 0.2 0.4 0.6 0.8 1 

ct/N 



M 0.2 



0, 





• 


1 §0W 

1^- — T , l 


0.2 0.4 


0.6 0.8 1 


1 , 

" — CLV 1 
- --■ OLV -J 


' ' '(d)- 

1 — =i , 



,1 I 1 , — - 1 --^-, -I— , .J 

0.2 0.4 0.6 0.8 1 

a/2N 



FIG. 5: Similar to Fig. 0] but the local map is the sinusoidal map and e = 1.3. The signifi- 
cance of HLM is strongly reduced in CLVs of nonhyperbolic dissipative systems whereas the HLM 
significance is comparable for CLVs and OLVs in nonhyperbolic Hamiltonian systems. 

r from 0.2 to 0.4 leads to a dramatic reduction of the significance of HLMs, in both CLVs 
and OLVs, for either dissipative system or Hamiltonian system. In the dissipative system 
the reduction of CLV significance as increasing r is much faster than the reduction of OLV 
significance, which leads to an increasing discrepancy between them. In contrast for the 
Hamiltonian system the significance of CLVs and OLVs is always comparable. Note, how- 
ever, that the tuning of r has no influence on the symmetry features of Lyapunov vectors. 



Owing to the different symmetry properties of CLVs and OLVs of Hamiltonian system, 
the largest S(k max ) is observed at different a values. This leads to a rather large difference in 
the static LV structure factors corresponding to the smallest positive Lyapunov exponents. 
We show in Fig. [8]two cases with different coupling strength e. As can be seen from the figure 
the CLV structure factor S(k) diverges quickly as k goes to zero while the OLV structure 
factor increases relatively slowly and even seems to saturate to a constant. Note that the 



10 




0.2 0.4 0.6 Oi 
OC/N 



FIG. 6: Influence of the hyperbolicity variation on the significance measure S(k max ) of HLMs in the 
dissipative system. The parameter r is (a) 0.2, (b) 0.3 and (c) 0.4 respectively, which corresponds 
to a decreasing hyperbolicity. The local map is the skewed tent map and e = 1.3. 
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FIG. 7: Influence of the hyperbolicity variation on the significance measure S(k max ) of HLMs in the 
Hamiltonian system. The parameter r is (a) 0.2, (b) 0.3 and (c) 0.4 respectively, which corresponds 
to a decreasing hyperbolicity. The local map is the skewed tent map and e = 1.3. 
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FIG. 8: Static LV structure factor corresponding to the smallest positive Lyapunov exponent of 
the Hamiltonian case. The local map is the sinusoidal map. The large difference between these 
specific CLVs and OLVs is due to the asymmetric locating of the largest value of S(k max ) of OLVs 
as shown in Fig. [5b. 
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FIG. 9: Variation of the significance measure S(k max ) of HLMs with the system size N for the 
Hamiltonian case with the sinusoidal map. The peak moves gradually to the spectral center point 
a/2N = 0.5 as increasing N. 



same parameter e = 0.6 was used in Ref. 14J (see Fig. 3 therein). 

Moreover, as increasing the system size N the asymmetrically located peak of S(k max ) 
for OLVs shifts towards the spectral center point a/2N = 0.5 as shown in Fig. |9l which 
indicates a gradual reduction of the discrepancy between CLVs and OLVs as approaching 
the thermodynamic limit. 



VII. CONJUGATE PAIR RELATION IN HAMILTONIAN SYSTEM 

By definition, CLVs are in general not mutually orthogonal as OLVs are. This difference 
has some interesting consequences in Hamiltonian systems. As reported in Ref. 
conjugate pair of OLVs with \^ = —\^ 2L ~ l ~ a ) has the symmetry that Sw a ' = ±5t> ( 2L_1 ~ a ) 
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and 5v^ = =F<W 2i 1 a \ Here 5u and 5v denote the coordinate and momentum parts 

of LVs respectively. The physical origin of this symmetry lies in the symplectic structure 

— T 

of Hamiltonian system. OLVs as the eigenvectors of the matrix H = lim^ +00 [M(t, 0) • 



M(t, 0)]5t are thus forced to have the observed symmetry. Examples of conjugate pairs of 
OLVs are shown in Fig. [TUJ 

As can be seen from the same plot, CLVs behave differently. The relation 5u^ = 

seems to work well instead. Note also that for CLVs 
the amplitude of the wave structure in the momentum part is much smaller than in the 
corresponding coordinate part. Such difference is also reflected in the profiles of S(k max ) in 
Fig. [TTJ For OLVs S(k max ) of the coordinate part and momentum part are mutual mirror 
images with respect to the spectral center a/2N = 0.5. For CLVs S(k max ) from either the 
coordinate part or the momentum part is roughly symmetric with respect to the center by 
itself. 

As going to the nonhyperbolic cases with the sinusoidal map as the local dynamics, the 
mentioned simple relations between Su and Sv of instantaneous Lyapunov vectors valid no 
longer. However, as can be seen from Fig. [T2J the conjugate pairs of CLVs have nearly 
identical k max and S^ a '(k max ), i.e. they are statistically indistinguishable. This interest- 
ing feature of CLVs is believed coming from the micro-reversibility of Hamiltonian system. 
Micro-reversibility means that for each trajectory from (u(0),v(0)) to (u(T),v(T)) there 
exists a reverse-time trajectory from (u(T),—v(T)) to (w(0), — v (0)). Under time reversal 
Lyapunov exponents change their sign and the conjugate pair of CLVs exchange their role 
for characterizing the stable and unstable directions. Owing to the ergodicity of Hamilto- 



nian system 



181 ] S^ a \k max ) for the pair of initial conditions (u(0),v(0)) and (u(T), —v(T)) 



are indistinguishable and this leads to the observed symmetry feature of CLVs. 

Besides the mentioned symmetry resulting from the general Hamiltonian property, the 
conjugate pair of CLVs in systems with continuous symmetries such as Eq. (j3J) have some 
unexpected interesting features. The angle 9 between a pair of CLVs is used to characterize 
their relation, with cos(#) = |e( Q ) • e^]. The contour plot of the quantity (cos(#)) is shown 
in Fig. [131 were (• • • ) means an average over time. 

It shows that for the highly hyperbolic cases, for instance Eq. (jSJ) with the special skewed 
tent map as the local dynamics, CLVs corresponding to near-zero Lyapunov exponents are 
nearly orthogonal to each other as expected. The fact is more evident in Fig. [T3J The near 
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FIG. 10: Instantaneous profiles of two conjugate pairs of CLVs and OLVs for the Hamiltonian case. 
The local map is the skewed tent map with e = 1.3 and r = 0.15. The system size used is N = 128. 
Note that OLVs and CLVs have different symmetry relations for the conjugate pair. 



orthogonal nature of those CLVs explains the observed similarity between CLVs and OLVs 
in Fig. [2] and H] for the current parameter setting. In contrast the conjugate pair of CLVs 
tend to the same orientation as approaching the zero Lyapunov exponents. With changing 
the local dynamics to the sinusoidal map the near orthogonal regime disappears completely 
whereas the qualitative behavior of the angle between conjugate pairs is hardly influenced. 



VIII. CONCLUSION AND DISCUSSION 



We have explored, by using simple models of coupled map lattices, the similarity and 
difference between CLVs and OLVs, especially with respect to hydrodynamic Lyapunov 
modes. For both Hamiltonian and dissipative cases, two different local maps were used to 
represent the typical situations with different degree of hyperbolicity. The dynamics of the 
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FIG. 11: Comparison of the dominant W3,V6 number k max and S( a \k 

max) obtained from the co- 
ordinate and momentum parts of OLVs (a,b) and CLVs (c,d) respectively. The local map of the 
studied Hamiltonian system is the skewed tent map with e = 1.3 and r = 0.15. 

case with the special skewed tent map is highly hyperbolic while the one with sinusoidal 
map is nonhyperbolic as most systems. In some sense the Hamiltonian system with the two 
local maps are corresponding to the often used hard-core system and soft potential system 
respectively. With the replacement of OLVs by CLVs the formerly detected long wave- 
length structure in Lyapunov vectors can be seen as well. Moreover the CLV X-k dispersion 
relation is linear for Hamiltonian system while quadratic for dissipative system as found for 
OLVs. The significance of HLMs as measured by the static LV structure factor changes 
differently for Hamiltonian and dissipative systems with the replacement of OLVs by CLVs. 
For Hamiltonian systems the significance of HLMs is always comparable for CLVs and OLVs 
independent of the variation of hyperbolicity as changing the local maps, besides that the 
OLVs with the most significant wave structure lie slightly away from the spectral center. 
Increasing system size tends to shift them back to the center. For dissipative systems the 
significance of HLMs is almost the same for CLVs and OLVs if the special skewed tent map is 
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FIG. 12: Similar to Fig. [11] but the local map is the sinusoidal map and e = 1.3. 

used as local map. Departing from such a highly hyperbolic situation the HLM significance 
of CLVs reduces much faster than that of OLVs. 

In the past there were already discussions regarding the symmetry of the conjugate pair of 
Lyapunov vectors in Hamiltonian system. It was found that owing to the symplectic feature 
of Hamiltonian systems the coordinate and momentum parts exchange their position for 
a conjugate pair of OLVs. A different symmetry is observed, however, for CLVs, namely 
that two CLVs in one conjugate pair are statistically indistinguishable. As discussed the 
physical origin of this seemingly unreasonable property is the microscopic reversibility, a 
general feature of Hamiltonian systems. For the specific issue of HLMs, it implies that the 
variation of HLM significance for CLVs is symmetric with respect to the spectral center. 
Besides that we found for CLVs that the HLM significance is much lower in the momentum 
part than in the coordinate part. 

For the highly hyperbolic cases with the special skewed tent map as the local dynam- 
ics CLVs behave very similar to OLVs as demonstrated by the position and height of the 
dominant peak of static LV structure factors. A direct monitoring of the mutual angle 
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FIG. 13: (a) Contour plot of (cos(0)) of CLVs for the Hamiltonian case. The local map is the 
skewed tent map with e = 1.3 and r = 0.15. Panel (b) shows the enlargement of the central part of 
(a). Note that in the regime A ~ CLVs are nearly mutual orthogonal besides that the conjugate 
pair of CLVs have a very small angle. 



between CLVs shows that for those corresponding to near-zero Lyapunov exponents the mu- 
tual angles are large and close to tt/2. An unexpected observation is that the angle between 
conjugate pair decreases to zero as approaching the spectral center. Such a feature persists 
as weakening the hyperbolicity. 

It is known that a dynamical system has two sets of OLVs, backward and forward ones. 
Only the backward OLVs, which can be calculated numerically via the standard method, 
are discussed in the main text. Similar results are expected for the forward OLVs except 
that they bear the similarity to a different part of CLVs compared to the backward OLVs. 
A related discussion can be seen in the appendix 

To mention that a comparison of CLVs and OLVs in systems with hard-core interactions 
is performed by Posch et al 2l|, which and the current contribution form a complementary 
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FIG. 14: The quantity (cos(0)) of angles between the conjugate pair (square) and neighbouring 
CLVs (circle) of the Hamiltonian system. The local map is (a) the skewed tent map with r = 0.15 
and (b) the sinusoidal map. The coupling strength is e = 1.3. Switching to the nonhyperbolic case 
the orthogonality between neighbours CLVs in the regime A ~ is broken. 



view of the topic to each other. 
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Appendix A: asymptotic and finite time Lyapunov exponents 

As can be seen from the definition and calculation algorithm in Sec. [TT] as well as from 
other sections CLVs and OLVs are different in many respects. In this appendix we would like 
to point out that the (asymptotic) Lyapunov exponents corresponding to CLVs and OLVs 
are identical but finite-time Lyapunov exponents (FTLEs) corresponding to these two sets 
of Lyapunov vectors are different in general. 

From the calculation algorithm we know that a fc-dimensional vector spanned k arbitrary 
offset vectors will approach asymptotically the most unstable /c-dimensional subspace, which 
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can be spanned by k Lyapunov vectors, either CLVs or OLVs, associated with the first k 
largest Lyapunov exponents. The growth rate of the volume V& of this k- dimensional 
subspace can be written as 

MUM) = i>fo(M a ) + r 1 T ln ' nS COS ^' 2 ? l (AD 

where A^^ti,^) is the growth rate of offset vectors along the z-th CLV, i.e. the z-th FTLE 
corresponding to this CLV and 9 % {t) is the angle between the i-th CLV and the subspace 
spanned by CLVs with index from to i — 1. Taking into account the mutual orthogonal 
nature of OLVs, the growth rate S fc can be expressed as well by using characteristics of 
OLVs as 

k-l 

S fc (ti,t 2 ) = ^A l (0) (t 1 ,t 2 ) (A2) 

i=0 

where X l , Jti,t2) is the growth rate of offset vectors along the i-th OLV, i.e. the z-th FTLE 
corresponding to this OLV. 

Combining Eq.( ]Al|) and (1A2j) yields a simple relation between FTLEs 

Aj 0) (t lj t 2 ) = Aj 0) (ti,t 2 ) (A3a) 

X {C) (tiM) = X\ 0) (tu+t2) + f^M^^I for z G [1, V- 1] (A3b) 

with N the dimension of the considered system. Since CLVs are in general not mutually 
orthogonal it is obvious from Eq. flA3bl) that FTLEs Av^ -j (t l5 1 2 ) and A/^ti,^) with i > 1 
are normally different. 

As approaching the limit t% — t\ = +oo the contribution of the second term in r.h.s. of 
Eq. (1A3bl) becomes negligible since the value of cos 9 is bounded. This implies 

Afofa.+oo) = X\ 0) (t u +oo) for any ie[0,JV-l], (A4) 

i.e. asymptotic Lyapunov exponents corresponding to CLVs and OLVs are identical. 
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Appendix B: transformation properties of Lyapunov exponents, CLVs and hyper- 
bolicity 

1. Invariance of Lyapunov exponents and covariance of CLVs 

We consider a dynamical system which is written as 

x = F(x) or x t+ i = F(xJ. (Bl) 

The time evolution of its trajectory can be expressed as 

x(t 2 ) =0(ti,t 2 )x(ti). (B2) 

Correspondingly the evolution of an infinitesimal perturbation vector with respect to the 
reference trajectory x(t) can be written as 

6x.{t 2 ) = M(t 1 ,t 2 )6x{t 1 ) (B3) 

with M = <90/<9x. Under a variable transformation T : x i— > y with 

y = T(x), (B4) 

the governing equation of infinitesimal perturbations becomes 

Sy(t 2 ) = M%,t 2 )Sy{t 1 ). (B5) 

Here the two variables 5y and <5x are related via a linear transformation L : <5x i— > 5y with 

5y = L<5x (B6) 

It is known that the linear transformation L is determined by the transformation T via 

L = D x T(x), (B7) 

where (D x T)y = dTi/dx-j. By using Eq.flEJ), flB3]) and (JESJ one can show that 

ML = LM. (B8) 

which means that M' and M are related via a similarity transformation M' = LML" 1 if L 
is invertible. 
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If e(ti) = e(x(ti)) is a CLV in the x-coordinate system, it satisfies the condition 



M(t 1 ,t 2 )e(t 1 ) = <7(ti,t 2 )e(t 2 ) (B9) 

with A = lim f2 _ t -3— In |cr(ti,t 2 )| being the Lyapunov exponent corresponding to this 
CLV. Multiplying the both sides of Eq. (1B9j) with L and using Eq. (1B8j) results in 

M'(ti, t 2 )Le(ti) = <t(£i, £ 2 )Le(£ 2 ). (BIO) 



Denoting e'(t) = Le(t)/||Le(t)|| one can reformulate Eq. flBlOl) as 



M'(ti, t 2 )e'(ti) = a(t 1; t 2 ) |^ e ^|| e'(t 2 ) = o>(ti,h)e%). (Bll) 

||L,e(tiJ|| 



Under the condition that 



lim —M l P¥4) = (B12) 

one can easily obtain that 

lim ___ ]n !</(*!, t 2 )| = lim ——ln^ti,^)!, (B13) 

t 2 — *l— >oo J 2 — Cl t2-*i-5>oo t 2 — jf^ 

which implies that (i) the unit vector e'(t) is a CLV in the y-coordinate system and it is 
related to e(t) via e'(t) = Le(i)/||Le(t)||; (ii) the asymptotic Lyapunov exponent associated 
with e'(t) is identical to the asymptotic Lyapunov exponent corresponding to e(t); (hi) the 
finite-time Lyapunov exponent X'(ti, t 2 ) = ^^-(ln \a(ti, t 2 )\ +ln ) in the y-coordinate 

system is different from the one \(ti,t 2 ) = ^^-ln |er(£i,i 2 )| in the x-coordinate system. 
For an invertible transformation T with the assumption that the reference trajectory is 



bounded in phase space one can easily show the boundedness of L [19[, i.e. 



L _ ||e|| < || Le || < L+||e|| (B14) 

for two constant L~ < L + < 00, which implies the validness of the condition stated in 
E, m . As diseusseLnQ these s on T can be weakened such that for non- 

invertible transformations T one can still get the invariance of Lyapunov exponents and the 
covariant transformation of CLVs. This is also confirmed by our numerical example below. 
As discussed already in Ref.[7|, via the transformation 

x\ = u[+ l - u\ (B15) 
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FIG. 15: Instantaneous profiles of CLVs of two systems related via variable transformation 
Eq. (|B15j ). As shown in panel (d) CLVs of the two systems are related via the relation given 
in Eq.(jB17]). 
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FIG. 16: Difference between finite-time Lyapunov exponents of two systems related via variable 
transformation Eq. (IB15j) . 

Eq.flB]) can be mapped to the following diffusively coupled CMLs 

x\ +1 = x\ + e[f{x\ +l ) + f^f 1 ) - 2/04)]. (B16) 
According to our above arguments CLVs of the two system are related via the transformation 

6xi a)l = c t {8uf - 8uf )l ) (B17) 



where c t is a time- dependent normalization factor. Numerical results shown in Fig. 
and [16] for a case with the sinusoidal map f(z) confirms our conclusions. Note that the 
transformation Eq.( lB15j) is non-invertible. 



2. Invariance of hyperbolicity under diffeomorphisms 



In viewing that e'(t) = Le(t)/||Le(t)|| one would expect that the absolute value of angles 
between CLVs is not invariant under the variable transformation. Whether the angle is zero 
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or not, i.e the feature of hyperbolicity, is expected to be preserved under diffeomorphisms. 
This conjecture is supported by the following arguments. 

If the variable transformation T given in Eq. (IB4p is a diffeomorphism, the corresponding 
transformation L of the perturbation in Eq. flB4j) would be an invertible linear transformation. 

Consider two CLVs ei(t) and e 2 (t) in the x-coordinate system and denote the correspond- 
ing CLVs in the ^-coordinate system as e[(t) and e' 2 (t). Since an affine transformation like L 
preserves the collinearity of points the angle between the transformed CLVs e[(t) and e' 2 (t) is 
zero if the angle between original CLVs ei(t) and e 2 (t) is zero, i.e. Z(ei(t), e 2 (t)) = implies 
/(e'^t), e 2 (t)) = 0. Similar arguments for the inverse L _1 leads to that /(e'^t), e' 2 (t)) = 
implies Z(ei(t), e 2 (t)) = 0. These properties indicate the preservation of the collinearity of 
CLVs under diffeomorphisms. 

Consider now two subspaces S\ and S2 spanned by two sets of different CLVs {e^} and 
{e^ 2 }. If the angle between the two subspaces is zero in one coordinate system it means 
that the two sets of CLVs are linearly dependent, i.e. Yli c « e Si + ^ e 5 2 = ^ ^ or cer tain 
constants c, and di. Preservation of collinarity under affine transformations implies that 
the corresponding transformed CLVs are linearly dependent as well, i.e. the angle between 
subspaces in the transformed coordinate system is also zero. Similarly one can show that if 
the angle between two subspaces is nonzero in one coordinate system it would be nonzero in 
other transformed coordinate systems, too. These arguments show that whether the angle 
between subspaces is zero or not is invariant under dffeomorphisms, i.e. the property of 
hyperbolicity is preserved. 

Similar to the discussion about the transformation properties of CLVs one can weaken 
the requirements on the transform T but rather the same conclusion about the hyperbolicity 
can be reached. A known example is the Kuromato-Sivashinsky equation. It can be written 
in two different forms as 

u t = u xx + u xxxx + u 2 x /2 (B18) 

or 

v t = v xx + v xxxx + vv x (B19) 

which are related via a non- invertible transformation v — u x . Numerical simulations show 
that the two forms have the identical hyperbolicity and details will be shown elsewhere [20 ] . 
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Appendix C: analytical calculation of CLVs and OLVs of a Hamiltonian system 



For the Hamiltonian system Eq.flSJ with the limiting case r = of the skewed tent map 
Eq.(J7J) one can calculate the CLVs and OLVs analytically. Consistence with numerical results 
presented in the main part of the paper can thus be checked. 

1. CLVs 

For the case r = in Eqs. (j5J7|) the time evolution of the infinitesimal perturbations is 
governed by 

Sf t+1= ( h+ ^ h ). S f t (CI) 

where 5Y t = {5u\ , 5uf, ■ ■ ■ , bu\\ 5vj, Svf, ■ ■ ■ , Sv^} is the offset vector in the tangent space 
and II, Dl denote the [L x L)-unit matrix and the discrete Laplacian, respectively. Notice 
that the fundamental matrix 

I I L + eD L I l \ 
M2 = f (C2) 
\ eD L I L J 

is time independent and thus the eigenvectors of M 2 are CLVs of this system. 

By using the eigenvectors e> a > of the matrix Dl the eigenvectors of the fundamental 
matrix M 2 can be constructed as {e^; c(k)(^ a ^}. The associated eigenvalues are 

y{k) + 2±y/rP(k)+4r,(k) 
fi±(k) = (C3) 

where T](k) = — 2e(l — cos/c) and the corresponding c(k) can be calculated as 

c±(k)=n±(k)-l-ri{k). (C4) 

The following properties of these eigenvectors/CLVs can be obtained: 

i) The corresponding eigenvalues satisfy fjL+(k)fjL-(k) = 1 which indicates the conjugate 
pair property A + = — A_ of Lyapunov exponents since A = In 

ii) The group of CLVs {e^ Q ^; c + (k)e^ a ^} corresponding to positive Lyapunov exponents are 
in general not orthogonal to CLVs {e^; C-(k)e^ a ^} corresponding to the negative branch of 
the Lyapunov spectrum, although members of either group are mutually orthogonal. This 
indicates that these eigenvectors/CLVs are not OLVs of that system. 
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iii) As approaching the spectral center A = 0, one has k — > 0. Two CLVs in a conjugate 
pair tend to be collinear, i.e. cos 6 — > 1, where 9 denotes the angle between that pair of 
CLVs. More precisely, as k — > 0, one has 

r)(k) « ek 2 , (C5) 

fi±(k) ^l±^k (C6) 

and 

c ± {k) « ±v^fc. (C7) 
Thus for a conjugate pair of CLVs {e^; \feke^°^ } and {e^; — A/e^e 1 ^} one has 

cos^^l-eA; 2 . (C8) 

In consistence to this, we reported in Sec. IVIII that for the cases with r close to as shown 
in Fig. [13] and [H] neighbouring CLVs are nearly orthogonal while those in a conjugate pair 
tend to be collinear as approaching the spectral center A = 0. 



2. OLVs 



Now we start to calculate OLVs of this system, which are eigenvectors of the matrix 
M^(Mj)" as n goes to infinity, where is the transpose of M 2 . 

The matrix Mi has the similar transformation M 2 = QAQ^ 1 , where the column vectors of 
Q are eigenvectors of M 2 and entries of the diagonal matrix A are corresponding eigenvalues 
H±(k) mentioned above. The matrix M 2 n (M 2 r ) n can thus be written as 

M 2 n (M 2 T )" = QA n Q~ 1 (Q~ 1 ) T A n Q T . (C9) 

Considering the orthogonal nature of , the discussion of eigenvalue and eigenvectors can 
be simplified by using the submatrix Q(k) and S(k) = M 2 ri (M 2 r )"(A;) related to the vectors 
&- a >{k). They are 




Q(k) =| I , (CIO) 

and 

S(k) = (c + (k)-c^k))- 2 \ a{k) m | . (Cll) 

\d(k) b{k) 

25 



4). 



where a = /4™(l + d) + 2/i>™ +/r?"(l + 4), b = 4/4" (l + <£ ) + 2c+c_ / u>™ + <£/£"(! + 
and d = c + /i+ n (l + c 2 _) + (c+ + c_)//"^ + c_/ri n (l + c+). The eigenvalues of the matrix S'(fc) 
can be obtained as 

wt) = "_±±±V£±SEME^. ( ci2) 

As the eigenvectors of the matrix M^(Mj) n can be constructed as {e^; pe^}, one can 
easily get 

p(jfe) = ^ = (fe-o)/2 + V(a + &) 2 -4(a6-^ 

Since /J>+(k) > 1 > fi-(k), as n goes to infinity one has — >• c + (/c) with the corresponding 
A + (fc) = ~ln|£ + (A;)| — > ln|/x(A;)|, which indicates that the OLVs associated with positive 
Lyapunov exponents are the same as the corresponding CLVs. 

In consistence to this, as shown in Fig. HU and [6^,, for cases with r close to OLVs and 
CLVs associated with positive Lyapunov exponents are very similar. 

A dynamical system has actually two sets of OLVs, namely backward and forward OLVs, 
which are eig envectors of the matrix M 2 n (M 2 T ) ri and (M 2 T ) n M 2 " as n goes to infinity, respec- 
tively. In above discussions the backward OLVs are used since they are the ones numerically 
obtained from the standard method 11]. For the forward OLVs one can do the similar 



calculations and the conclusion is that a half of the forward OLVs are the same as the CLVs 
associated with negative Lyapunov exponents. 



[1] J. P. Dorfman, An Introduction to Chaos in Nonequilibrium Statistical Mechanics (Cambridge 
University Press, Cambridge, 1999); P. Gaspard, Chaos, Scattering, and Statistical Mechanics 
(Cambridge University Press, Cambridge, 1998). 

[2] E. Kalnay, Atmospheric Modeling, Data Assimilation and Predictability (Cambridge Univer- 
sity Press, Cambridge, 2002). 

[3] H.A. Posch and R. Hirschl, in Hard Ball Systems and the Lorentz Gas, ed. D. Szasz, (Springer, 
Berlin, 2000), p. 279. 

[4] J.-P. Eckmann and O. Gat, J. Stat. Phys 98, 775 (2000); A.S. de Wijn and H. van Beijeren, 
Phys. Rev. E 70, 016207 (2004); T. Taniguchi and G.P. Morriss, Phys. Rev. Lett. 94, 154101 
(2005). 

[5] S. McNamara and M. Mareschal, Phys. Rev. E 64, 051103 (2001). 

26 



[6] H.L. Yang and G. Radons, Phys. Rev. E 71, 036211 (2005); G. Radons and H. L. Yang, 

|arXiv:nlin.CD/0404028| 
[7] H.L. Yang and G. Radons, Phys. Rev. E 73, 016202 (2006). 
[8] H.L. Yang and G. Radons, Phys. Rev. Lett. 96, 074101 (2006). 
[9] A. Pikovsky and A. Politi, Nonlinearity 11, 1049 (1998). 
[10] V.I. Oseledec, Trans. Moscow Math. Soc. 19, 197 (1968). 

[11] G. Benettin, L. Galgani and J. M. Strelcyn, Phys. Rev. A 14, 2338 (1976); I. Shimada and 

T. Nagashima, Prog. Theor. Phys. 61, 1605 (1979). 
[12] I. Goldhirsch, PL. Sulem and S.A. Orszag, Physica 27D, 331 (1987). 
[13] S.V. Ershov and A.B. Potapov, Physica D 118, 167 (1998). 
[14] F. Ginelli, et al., Phys. Rev. Lett. 99, 130601 (2007). 

[15] I.G. Szendro, D. Pazo, M. A. Rodriguez, and J. M. Lopez, Phys. Rev. E 76, 025202(R) (2007). 

[16] K. Kaneko, Prog. Theor. Phys. 72, 480 (1984). 

[17] H. L. Yang and G. Radons, Phys. Rev. Lett. 100, 024101 (2008). 

[18] For those Hamiltonian systems which are nonergodic one can always find a decomposition of 
the invariant measure to ergodic components and we expect that the two trajectories talked 
are in the same ergodic component. 

[19] R. Eichhorn, S.J. Linz and P. Hanggi, Chaos, Solitons and Fractals 12, 1377 (2001). 

[20] H.L. Yang and G. Radons, in preparation. 

[21] H.A. Posch, et al, a talk given in the international workshop "Exploring complex dynamics 
in high-dimensional chaotic systems: from weather forecasting to oceanic flows", Dresden 
Germany, 25-29 Jan. 2010, organized by J.M. Lopez, A. Pikovsky and A. Politi. 



27 



